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Abstract. In this paper wc study configuration spaces of hard spheres in a 
bounded region. We develop a general Morse-theoretic framework and show 
that mechanically balanced configurations play the role of critical points. As 
an application, we find the precise threshold radius for a configuration space 
to be homotopy equivalent to the configuration space of points. 

1. Introduction 

Configuration spaces of n points in R'' are well studied [4]. In this article we are 
interested in a natural generalization, configuration spaces of non-overlapping balls 
in a bounded region in M'^. 

Besides their intrinsic mathematical interest, the study of these spaces is mo- 
tivated by physical considerations. For example, in statistical mechanics "hard 
spheres" (or in two dimensions "hard disks") are among the most well-studied 
modcils of matter. Computer simulations suggest a solid-liquid phase transition for 
hard spheres [17], but this is not well understood mathematically. 

A number of papers in statistical mechanics have explored the hypothesis that 
underpinning phase transitions are changes in the topology of the underlying con- 
figuration space or equipotential stibmanifolds [24, 16, 1, 9]. Franzosi, Pettini, and 
Spinelli show that under fairly general conditions (smooth, finite-range, confining 
potentials), the Helmholtz free energy cannot pass through a phase transition unless 
there is a change in the topology of the underlying configuration space [11, 10]. This 
theorem unfortunately does not apply to configuration spaces of hard spheres, since 
the potential function is not smooth — but the Morse-theoretic methods developed 
here may be a step in the direction of extending it to include hard spheres. 

Several other papers have investigated configuration spaces as models of motion 
planning for robots [8, 13]. For example, Farber's "topological complexity" can be 
thought of as measuring the difiiculty of designing an algorithm for navigating the 
space. As Deeley recently pointed out when he studied "thick particles" on metric 
graphs, the assumption that robots are points is not physically realistic, and giving 
the points thickness wildly complicates the topology of the underlying configuration 
space [6]. 

Let Bhc a bounded region in M.'^. Define Conf (n, r) to be the configuration space 
of n non-overlapping balls of radius r in B. We are especially interested here in 
understanding when the topology changes if n is fixed and r is varying . First we 
consider the extreme cases. For r sufficiently small, one expects that Conf (n, r) is 
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homotopy equivalent to the configuration space Conf (n) of n distinct points in B 
— for a survey of configuration spaces of points see Colien [4] . On the other hand 
for r sufficiently large, Conf (n,r) is empty. Indeed finding the smallest such r is 
the sphere packing problem in a bounded region — see for example Graham et al. 
[15, 2, 18, 19] and Melissen [22, 21]. 

In this note we develop a Morse-theoretic framework which provides a necessary 
condition for the topology to change — mechanical balanced configurations play 
the role of critical points (and submanifolds) . As an illustration of the method, we 
find the precise threshold radius below which Conf (n, r) is homotopy equivalent to 
Conf (n). 

2. Tautological Morse function 

Fix n, and define Conf (n) to be the set of ordered n-tuples of distinct points in 
a bounded domain B CR"^: 

Conf (n) = {x = {xi,. . . ,Xn) \ Xi € B,Xi Xj for i ^ j}. 

As an open subset of M''", Conf (n) has the structure of a smooth manifold. Let 
r : Conf (n) ^- M be defined by 

(1) t{x) := min ( ^ min(i(a;i,a;j),min min ) , 

\^ 2 i^o i p^dB J 

where dB denotes the boundary of B. We call r the tautological function. Then 
by definition the configuration space of n balls of radius r in S is given by 

Conf (n, r) ~ [r, oo). 

This observation suggests using a "Morse" -type theory of t to study the topology 
of Conf(n,r) and especially how the topology changes as r varies. One obvious 
trouble on that route is the fact that r typically is not smooth, so that we need a 
general framework which allows us to work with non-smooth functions. 

In the next section we will discuss the properties of min-type functions, that is 
functions on a manifold M given as the minimum of a parametric family of real 
valued functions 

t{x) := min/(p, x),x £ M,p € P, 

V 

where P is a compact parameter space, and / is continuously differentiable in x 
for every fixed p G P. We note that the function given by (1) falls within this 
category, if one considers as P the disjoint union of the discrete set corresponding 
to pairs («, j), 1 < i < j < n and of n copies of the boundary, formed by the pairs 
ihp), I < i < n,p G dB. 

It should be remarked that the Morse-type theory of the mAn (or even min-max) 
type functions has appeared in the literature (compare [20, 3, 12]), but in a much 
more restrictive context (with essentially finite parameter space). 

3. MiN-TYPE Morse theory 

Let us start with some notation and definitions. 

For a manifold M and function / : M — >■ K, let M'^ denote the superlevel set at 
c, i.e. = /~^[c, oo). 
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We say that a function h : [s, t) — ?• M, is increasing with speed at least t; > if 
h{t') — h{s') > v{t' — s') for any s' < t' in the interval {s,t). We note that such a 
function does not have to be even continuous. 

We record for later use the following (immediate) result: 

Lemma 3.1. Let f : M x R ^ R be a continuous proper function, (strictly) 
increasing along each fiber {x} x M. Then the fiber-wise inverse <p : {x, c) i— >■ iid{t : 
f{x,t) = c) is continuous on M x [ci,c2] for any interval [ci,C2] which belongs to 
the ranges of all the functions fx{-) '■— f{x, •). 

Proof. As / is proper, then the t-projection of the preimage of a compact interval 
[ci,C2] is compact as well. Further, if a sequence (xi,Cj) converges to (a;*,c*), yet 
ti := (j){xi,Ci) fails to converge to f» := 0(a;,,c,), we can, using the boundedness 
of the sequence (ti), choose a subsequence such that along it tj t ^ t^.. By 
continuity, Yim f{xj,tj) — f{x^,t) (as Xj converge to cc*). The fact that f{xj,tj) = 
Cj implies that f{x^,t) equals c,, which by the continuity of / also equals f{x^,t^). 
This contradicts the assumption that / increases fiber-wise. □ 

For a smooth vector field y on M we will denote the time t shift along the 
trajectories of V as SY ■ We will say that the function / increases along the tra- 
jectories of V with non-zero speed, if for some common v > 0, and for all x € M, 
hx ■ 1 1-^ f{SYx) increases with speed at least v. 

Lemma 3.2. Let M be a smooth manifold and / : M — >• R a continuous function, 
such thai A'/" is compact. Suppose that M admits a, smooth vector field V non- 
vanishing on f~^[a,b], and such that f is increasing along the trajectories ofV on 
the set f~^[a,b] with non-zero speed. Then is a deformation retract of M'^. 

Lemma 3.2 can be seen as a generalization of Theorem 3.1 in [23] for non-smooth 
functions /. We remark here that one can drop here the non-zero speed condition, 
requiring only that / is increasing along the trajectories of V, but we do not need 
this strengthened form in this paper. 

Proof. For x G M" set a partially defined function on x M by g{x, t) := f{SYx). 
The non-zero speed condition implies that [a, 6] is in the range of g{x,-) for any 
X € /~^([a, b]), and together with the compactness of M" implies that g is proper. 
Hence, by Lemma 3.1, 4>{x, c) = inf (t : f{SYx) > c) is well defined and continuous, 
as well as (/)(x, c) = max(^(x, c), 0). As for any c < f{x), (j){x, c) =0, (f) vanishes on 
M'' X [a,b]. 

Now define the homotopy 

H -.M" X [0, 1] M" 

as 

H : {x,t) ^ S'^(^,(i_^)a+^6)a;. 

Continuity of (j) implies continuity of H: the facts that H{x, 0) = idM", H{x, t)\m'' = 
idj^jb for < T < 1 and H{x, 1) € are immediate. □ 

3.1. Regular values of min-type functions. Next we will use the fact that the 
tautological function r is the minimum of a compact family of smooth functions. 

We want to establish conditions when a vahie c is topologically regular, that is for 
which there exists some e > such that 71^^+^ is a deformation retract of M°~*. 
We give a general condition for topological regularity, as follows. 
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Let P be a compact metric space, M a compact smooth manifold with boundary 
and 

/ : P X M ^ M 

a continuous function such that the x-derivative of / (that is, the gradient of fp, 
where fp{x) = f{p,x)) is continuous on P x M. We will be talking of P as the 

parameter space. 

We denote by r := minpgp/j, the min- function of the family /. The set N C 
P X M defined by 

N:={{p,x):f{p,x)=T{x)} 

is compact, and the slices 

N,:^{peP:{p,x)eN} 

are upper semi-continuous: for any x € M and any open neighborhood U D 
there exists an open neighborhood Ux 3 x such that for x' e Ux, N^i C UN^- 

Next wc show that if one can perturb each x to increase r then we can do so 
globally with a minimum speed. 

Lemma 3.3. Assume that for any x G M , there exists a tangent vector Vx S T^M 
such that Lv^ fp>0 for all p € N^. Then 

• for some positive v there exists smooth vector field V on M such that 

Lvfp > V > in som,e open vicinity of N , and 

• along the trajectories of V , the mhi- function r increases with speed at least 

V. 

Proof. For any x G M, wc can extend the vector S T^M to a smooth vector 
field on M (which wc still denote as Vj-), such that Ly^f > in some open vicinity 
UNx X Ux of Nx X {a;}. By compactness, there exists a finite collection of points 
{xj} in M such that the open sets Ui := UN^^ x Uxi cover N (and the open sets Uxi 
cover M), and u > such that Lv^fp > v on Ui (here Vi := 14; )• Using a partition 
of unity we arrive at the first conclusion. The second conclusion is immediate. □ 

For X G M consider the intersection of the open half-spaces 

Hxip) ■■= {vGTxM ■.{dfp\x,v)>0} 

over all p £ N^. This is an open convex cone 

C° := fl Hxip) 

in TxM. 

Upper semicontinuity of implies lower semicontinuity of C°; for any x G M 
and any open set V C TM intersecting C°, there exists an open neighborhood 
Ux 9 X such that for x' G Ux, C°, intersects V. In particular, if C° is non-empty, 
it remains such in a vicinity of x. 

Combining Lemmata 3.2 and 3.3 we obtain the following 

Corollary 3.4. If the cones C° are non-empty over the level set r~^(c), then c is 
topologically regular. 

For general mm-typc functions this is essentially the best possible condition 
for the regularity of the critical values. If the functions fp are quasi-convex, i.e. 
have convex lower excursion sets {fp < c}, then Corollary 3.4 can be considerable 
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strengthened. Thus, one can show (we will do it in a follow-up paper) that a critical 
value is topologically regular, if for all points at the level set, the intersection of 
the dosed half-spaces is a cone over a contractible base. This observation relies on 
stratified Morse theory due to Goresky and MacPherson[14], but is in a nutshell 
close to the elementary result used by Connelly in his work on the existence of 
continuous "unlocking" deformations of hard ball configurations, see [5] . 

Corollary 3.4 implies that unless the level set of the tautological function T~^(r) 
contains a point x with C° = 0, the homotopy type of Conf (n, r) is locally constant 
at r. 

By Farkas' lemma, the emptiness of the cone C° implies that there exists a finite 
collection of points Pi G Nx,i = 1, . ■ . , I < dim M + 1, and positive weights w, > 
such that 



In our hard spheres setting, the vanishing of the convex combination (2) has a 
clear geometric interpretation. 

For X G Conf (n, r) , define a stress graph of a; to be a graph embedded in R'^ whose 
vertices are the points xi, . . . , x„ and boundary points y G dB where d{xi, y) = r 
for some ?'. The edges are the pairs {xi, Xj} where d{xi, Xj) = 2r and {xi, y} where 
d{xi,y) = r. Each edge k is assigned a positive weight Wk. The points Xi are 
referred to as internal points and the points y are referred to as boundary points. 
Wc interpret this graph as a system of mecha,nical stresses, with (repulsive) forces 
acting on the endpoints of a segment k equal to Wk times the unit vector in the 
direction of k. Call the mechanical stresses acting on boundary points boundary 
mechanical stresses. Call a connected component trivial if it consists of a single 
point. Call x trivial if r(a;) has no edges. 

A stress graph is said to be balanced if it satisfies the following condition. 

• The mechanical stresses at each internal point sum to zero.^ 

• The boundary mechanical stresses on each connected component sum to 
zero. 

Say that the configuration is x is balanced if it has a balanced stress graph. 

Call an internal point isolated if it is not in the boundary of any edges. For 
each point Xi call the intersection of the stress graph with the points on the sphere 
d{xi,x) = r kissing points of Xi. Call a kissing point that is also a boundary point 
a boundary kissing point. 

Lemma 4.1. Assume that x € Conf (n,r) is balanced. Then 

(1) each non-isolated internal point is in the convex hull of its kissing points, 
and 

(2) each non-trivial connected component is contained in the convex hull of its 

boundary kissing points. 

Now consider x & Conf (n) that is a critical point of r with critical value r. 
In (2), the parameters Pi correspond either to the pairs of touching hard spheres, 
d{xai,Xbi) = 2r, or to the hard sphere x^ touching the boundary, d{xci,yi) = r, 

^This result is similaj: to the necessary conditions for "locking" , see e.g. [5] . 



(2) 




4. Critical points and stress graphs 
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at a point yt e dB. Let T{x) be the corresponding stress graph for x with weights 
given by the coefficients Wi in (2). 

Theorem 4.2. If x & Conf (n) is a critical point of t with critical value r, then x 

is balanced and nontrivial as a point in Conf (n, r). 

Proof. Consider x G Conf (n,r) and consider T{x). From (2) it follows that the 
mechanical stresses at each internal point sum to zero. 

The sum of mechanical stresses in each connected component equals the sum of 
mechanical stresses on internal points and th(> sum of (^xtcrnal mechanical stresses. 
From (2) and the first observation it follows that the; boundary mechanical stresses 
on each connected component sum to zero. 

Finally, since the sum in (2) is nontrivial, x is nontrivial. □ 

5. Hard spheres in a box 

Consider now in more detail the case of hard spheres in a rectangular box with 
sides L := Li < . . . < L^i, given, for dcfinitivencss sake, by 

B = {0 < fm < L„i, m = 1, . . . ,d}. 

(Here {fm} is the orthonormal coordinate system on WK) 

5.1. Initial interval. We show that Theorem 4.2 implies a lower bound on the 
length of the initial interval of values of r, where the homotopy type of Conf (n, r) 
remains constant. 

Theorem 5.1. For the rectangular box B, there are no critical values of r in 

(0,L/2n), and therefore, 

Conf (n,r) ~ Conf (n) 

for r < L/2n. 

Proof. Assume that x G Conf (n) is a critical value for r with critical value r. Then 
by Theorem 4.2, x is balanced and nontrivial. A connected component of T{x) 
contains at most n internal points and thus has diameter at most 2nr. Since x is 
nontrivial, it has at least one nontrivial connected component. It is contained in 
the convex hull of its boundary points, so it contains at least one boundary point. 
Since the boundary mechanical stresses of this connected component sum to zero, 
it must contain a pair of boundary points from opposing faces. Thus the diameter 
of this connected component is at least L. Therefore r > L/2n. □ 

We remark that the balanced stress graph of minimal diameter is not necessarily 
a segment, for non- rectangular boxes. For example, for the "concave triangle", it 
is a cone over three points, see Figure 1. 

5.2. First perestroika. A natural question is now to ask. whether there is a topol- 
ogy change as nr goes above the minimal length of the stress graph. We concentrate 
in the rest of the note on the case of the rectangular box with the shortest side of 
length L, and will investigate, whether 

i : Conf (n,r') Conf (n,r),r' = L/2n + e,r = L/2n- e 

is a homotopy equivalence, for small enoiigh e. 

We argue that it is not, by presenting explicit nontrivial {dn — n — (i)-cycles in 
ker(ffi) c Ha„-n-d{Conf{n,r'),Z). 
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Figure 1. Minimal length stress graph 

Indeed, let < e < L/2n{n — 1) (so that {n — 1) disks of radius r' would fit 
within the box when arranged in a vertical column, and n would not), and consider 
the set Se of n-point configurations in B given by the conditions 

• xi fixed is at distance r' from the center of the face {/i = 0}, 

• \xi-i-i — Xi\ = 2r' for i — 1, . . . n — 1, and 

• Xn is at distance r' from the face {/i = Li}. 

In other words, we consider the configurations for which the n disks touch each 
other and the opposite horizontal faces, forming a chain, see Figure 2. 




Figure 2. A configuration in S^. 

An immediate computation shows that is diffeomorphic to a (nd — n — d)- 
dimensional sphere. Orient it in some way, obtaining a class s e Hnd-n~d{Cojif(n, r')). 

Next we show that s is nontrivial by constructing a cohomology class with which 
it has a nontrivial pairing. 

Consider now the set S of configurations in S" given by 

• all points Xi, . . . , a;„ have the same coordinates fs, ■ ■ ■ , fd', 

• all points X2, ■ ■ ■ ,Xn have the same coordinates f2, ■ ■ ■ , fd] 
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• the fl coordinates of X2, ■ ■ ■ ,Xn satisfy 
hixi) > r;fi{xi+i) - fi{x,) > 2r,forz = 2, 



,n- l;fi{xn) < L-r, 



and 



• /2(a;i) < f2{x2)- 



In other words, the configuration consists of n — 1 verticaUy ahgned nonover- 
lapping r-disks constrained to have the same (/i, /2)-plane (with the same /a, . . -d 
coordinates as the disk xi), see Figure 3. 



The conditions above are given by a finite coUection of hnear equahties and 
inequahties, and therefore define a convex polyhedron of dimension d + n. The 
boundary of this polyhedron is in B" — Conf(n, r'), whence, upon orientation it 
defines a relative class a e Hn+d{l3^^ , B''^ — Conf (n,r')). 

We notice that the space B" of n-tuples of points in B can be embedded into 
the nd-dimensional sphere S'"'^ (consider a large ball containing B" and contract 
its boundary to a point). By excision and the long exact sequence for a pair, 



i?„+d(S",S" - Conf(n,/)) = i?„+rf(5"^5"'^ - Conf(n,r')) = i?„+rf-i(5"'^ - 



Conf (n, r')). By Alexander duality the class a can be identified with a class (which 
we stih denote by a) in i?"''-"-'^(Conf (n, r')). 

Lemma 5.2. The pairing between the classes s and a is non-trivial: s ■ a = ±1. 

Proof. Indeed, the manifolds and S intersect transversally at a single point. □ 

As one can observe, there exists a retraction of to a point staying within 
Conf (n,r), implying that the class s is in the kernel of Hi. Indeed, we first can 
reduce all the distances between by shrinking the differences between the adjacent 
chain centers so that 



and xi remains fixed (clearly, this homotopy keeps the configuration in Conf (n, r)). 
Then one can pull all the vectors (x^+i — Xi) so that they point vertically upwards^ 
(as not one was initially pointing downwards). 




Figure 3. A configuration in E. 



2^1+1 ~ 2;j I— > 



- (1 - t)r' 



{xi+i - Xi),t^ [0, = 1, . . . ,n - 1 



"^We think of fi as height. 
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For each permutation tt of indices 1, . . . ,n one obtains different classes Stt, and 
one can easily see that the pairing with the corresponding n\ classes o"^ is non- 
degenerate (because corresponding and S's are all geometrically distinct). Hence, 
the rank of the kernel of Hi is at least n!. 

We notice that the sphericity of the set 5^ does not depend on the fact that the 
stress graph is a chain. For the configuration on Figure 1, the corresponding set is 
diffeomorphic to a sphere as well: it is just the corollary of the first critical value 
coming from a topologically Morse critical point, compare [20]. 

5.3. Betti numbers. We can also compute how the Betti numbers change across 

the first threshold. Set = L/2n, and note that as the tautological function is 
semi-algebraic for semi-algebraic regions, its critical values are isolated. As the 
only balanced stress graphs in the case of a rectangular domain are the chains 
spanning the shortest dimension, for some small e there are no other critical values 
in (r, - e,n + e). 

It is well known [4] that the configuration space Conf (n) of n (labeled) points in 
M** has Poincare polynomial 

P{t) :=^/3,f 

i>0 
n-1 

= 11(1+^*'-^) 

i=l 

= ! + ■■■ + {n- l)!if„_it("-2)(''-i) + (n - l)!*^"-^)^'^-^), 

where 

71-1 

1=1 

This tells us the Betti numbers of Conf (n, r« — e) , since we have already shown 
that Conf (n,r* — e) is homotopy equivalent to Conf(n). We wish to compute the 
Betti numbers of Conf (n,r* -I- e). 

Let N = (n— l)(d— 1). As we shrink the disks across the critical value r, = L/2n, 
to the configuration space we attach kn\ cells of dimension N, where k is the 
largest number such that Lk = L, whose boundaries are representatives for the 
homology classes s defined in Section 5.2. Each of these cells either increments /Jat 
or decrements /3jv-i- 

The first observation is that 

^j[Conf (n, n + e)] = ft [Conf (n, r* - e)] 

for z < — 2. As (the proof to appear in a follow-up paper) one can show that 
/3i(Conf (n, r)) = for i > N and r > r*, so in particular 

ft[Conf(n,r*-Fe)] =0 

for i > N. Thus (n — 1)! of the A''-cells increase /3n. This leaves only y8jv-i to 
compute. 

Every A^-cell that does not contribute to /Jat decreases Pn-i- Since we know 
that kn\ cells are added, (n — 1)! of them contributing to /3jv, we have 

/3jv- 1 [Conf (n,r* -|- e)] = /3jv-i[Conf (n, - e)] + kn\ - {n - 1)1 , 
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and since 

/3jv-i[Conf (n, r* 

we have 

/3Ar_i[Conf (n,r* + e)] 

6. Concluding remarks 

In a future article wc will discuss non-dcgcneracy of critical points, which is 
closely related to the question of making our necessary condition for a change in 
the topology sufficient. We also discuss defining and computing the index of critical 
points, and especially investigate more of the asymptotic properties of Conf(n,r) 
as n — )■ oo. In particular we obtain bounds on the rate of growth of Betti numbers. 

An important special case for which little seems known is: What is the threshold 
radius r = r{n) for connectivity of Conf (n, r)? This is an important question phys- 
ically, since for example ergodicity of any Markov process hinges on connectivity 
of the state space. Diaconis, Lebeau, and Michel noted that r < c/n is sufficient 
to guarantee connectivity of Conf(n, r) [7] and this is best possible for certain re- 
gions. It would be interesting to know if connectivity of the configuration space 
ever extends into the thermodynamic limit, i.e. are there any bounding regions so 
that Conf (n, r) is connected for r < Cn~^^'^ and some constant C > 0? 

Acknowledgments 

We thank AIM and for hosting the workshop on, "Topological complexity of 
random sets" in August 2009, where we started discussing some of these problems. 
Y.B. was supported in part by the ONR grant 00014-11-1-0178. M.K. thanks IAS 
for hosting him this year and Robert MacPherson for several helpful conversations. 

References 

[1] Luca Angclani, Lapo Casctti, Marco Pcttiiii, Giancarlo Ruocco, and Francesco Zamponi. 
Topology and phase transitions; From an exactly solvable model to a relation between topol- 
ogy and thermodynamics. Phys. Rev. E, 71(3):036152, Mar 2005. 

[2] David W. Boll, Jerry Donovan, Ronald L. Graham, and Boris D. Lubachevsky. Improving 
dense packings of equal disks in a square. Electron. J. Combin., TiResearch Paper 46, 9 pp. 
(electronic), 2000. 

[3] L. N. Bryzgalova. The maximum functions of a family of functions that depend on parameters. 

Funktsiorial. Anal, i Prilozhen., 12(l):66-67, 1978. 
[4] Frederick R. Cohen. Introduction to configuration spaces and their applications. In Braids, 

volume 19 of Led. Notes Ser. Inst. Math. Sci. Natl. Univ. Singap., pages 183-261. World 

Sci. Publ., Hackcnsack, NJ, 2010. 
[5] Robert Connelly Rigidity of packings. European J. Combin., 29(8):1862 1871, 2008. 
[6] Kenneth Deeley. Configuration spaces of thick particles on a metric graph, submitted, 2011. 
[7] Persi Diaconis, Gilles Lebeau, and Laurent Michel. Geometric analysis for the Metropolis 

algorithm on Lipschitz domains. To appear in Invent. Math. 
[8] Michael Farber. Invitation to topological robotics. Zurich Lectures in Advanced Mathematics. 

European Mathematical Society (EMS), Zurich, 2008. 
[9] Michael Farber and Viktor Fromm. Homology of planar telescopic linkages. Algebr. Geom. 

TopoL, 10(2):1063-1087, 2010. 
[10] Roberto Franzosi and Marco Pettini. Topology and phase transitions II. Theorem on a nec- 
essary relation. Nuclear Physics B, 782(3) ;219 - 240, 2007. 
[11] Roberto Franzosi, Marco Pettini, and Lionel Spinelli. Topology and phase transitions I. Pre- 

Uminary results. Nuclear Physics B, 782(3):189 - 218, 2007. 



r in^l)\Hn-i ■.d = 2 
\Q ■.d>3, 

( {Hn_i + kn - l){n - ly. :d = 2 
1 {kn-l){n-iy. :d>3. 



MIN-TYPE MORSE THEORY FOR CONFIGURATION SPACES OF HARD SPHERES 11 



[12] V. Gcrshkovich and H. Rubinstein. Morse theory for Min-typc functions. Asian J. Math., 
1(4);696-715, 1997. 

[13] Robert Christ. Configuration spaces, braids, and robotics. In Braids, volume 19 of Led. Notes 
Ser. Inst. Math. Sci. Natl. Univ. Singap., pages 263-304. World Sci. Publ., Hackensack, NJ, 
2010. 

[14] Mark Goresky and Robert MacPherson. Stratified Morse theory, volume 14 of Ergebnisse 

der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)]. 

Springer- Verlag, Berlin, 1988. 
[15] R. L. Graiiam and B. D. Lubachevsky. Repeated patterns of dense packings of equal disks in 

a square. Electron. J. Combin., 3(l):Research Paper 16, approx. 17 pp. (electronic), 1996. 
[16] Paolo Grinza and Alcssandro Mossa. Topological origin of the phase transition in a model of 

DNA dcnaturation. Phys. Rev. Lett., 92(15);158102, Apr 2004. 
[17] Hartmut Lowen. Fun with hard spheres. In Statistical physics and spatial statistics (Wup- 

pertal, 1999), volume 554 of Lecture Notes in Phys., pages 295-331. Springer, Berlin, 2000. 
[18] B. D. Lubax;hevsky and R. L. Graliam. Curved hexagonal packings of equal disks in a circle. 

Discrete Comput. Geom., 18(2):179-194, 1997. 
[19] Boris D. Lubachevsky, Ron L. Graham, and Prank H. Stillinger. Patterns and structures 

in disk packings. Period. Math. Hungar., 34(1-2): 123-142, 1997. 3rd Geometry Festival: an 

International Conference on Packings, Coverings and Tilings (Budapest, 1996). 
[20] V. I. Matov. Topological classification of the germs of functions of the maximum and minimax 

of families of functions in general position. Uspekhi Mat. Nauk, 37(4(226)):167— 168, 1982. 
[21] Hans Melissen. Densest pa<;kings of eleven congruent circles in a circle. Geom. Dedicata, 

50(1): 15-25, 1994. 

[22] J. B. M. Melissen. Optimal packings of eleven equal circles in an equilateral triangle. Acta 

Math. Hungar., 65(4):389-393, 1994. 
[23] J. Milnor. Morse theory. Based on lecture notes by M. Spivak and R. Wells. Annals of 

Mathematics Studies, No. 51. Princeton University Press, Princeton, N.J., 1963. 
[24] Ana C. Ribeiro Teixeira and Daniel A. Stajriolo. Topological hypothesis on phase transitions: 

The simplest case. Phys. Rev. E, 70(1):016113, Jul 2004. 

Departments of Mathematics and ECE, UIUC, Urbana, IL 
E-mail address: ymbSuiuc.edu 

Department of Mathematics, Cleveland State University 
E-mail address: p.bubenik@csuohio.edu 

School of Mathematics, Institute for Advanced Study, Princeton NJ 08540 
E-mail address: mkalile@math.ias.edu 



